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I. INTRODUCTION 
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An analysis of the eigenstates of a symmetry-broken spin-boson Hamiltonian is performed by 
computing Bloch and Husimi projections. The eigenstate analysis is combined with the calculation 
of absorption bands of asymmetric dimer configurations constituted by monomers with nonidentical 
excitation energies and optical transition matrix elements. Absorption bands with regular and 
irregular fine structures are obtained and related to the transition from the coexistence to a mixing 
of adiabatic branches in the spectrum. It is shown that correlations between spin states allow for 
q^-. an interpolation between absorption bands for different optical asymmetries. 
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QN . The spectral properties of the spin-boson Hamiltonian have frequently been adressed in the past and in the last 
years investigations on the spectral features of classical nonintegrability in the quantum case were performed (see e.g. 
jj] Q and references therein). The interest in the spectral and dynamical properties of the spin-boson Hamiltonian 
is due to its numerous applications which include a large variety of phenomena in molecular and solid state physics 
with one realization being the optical and transfer properties of excitonic-vibronic coupled dimers. Thus the spectral 
and related optical properties of the corresponding dimer Hamiltonian have intensively been investigated over the 
t-H \ years, as examples we point to the work exposed in |^,^|. Most of these investigations were based on a Hamiltonian 
t-H , corresponding to a symmetric dimer configuration. However, from the point of view of a general asymmetric situation 
ON ' the symmetric case is singular, e.g. the quantum states possess parity properties in the symmetric case, which are not 
present in general. It should be noted that dimers often constitute subunits of more complicated molecular aggregates 
with an asymmetric structure, which breaks the symmetry of the dimer configuration (see e.g. 0). Therefore the aim 
Oh, of this paper is to go beyond the symmetric case and to connect the properties of the eigenstates of a generalized and 
^Ij ■ symmetry-broken spin-boson Hamiltonian with the line spectrum of absorption bands of asymmetric excitonic-vibronic 
coupled molecular dimers. 

We combine a phase space analysis of the eigenstates, which is based on the method of Husimi projections [||, 
with the calculation of the matrix elements of optical transitions. Husimi projections have repeatedly been used 
in order to establish the features of nonintegrability of classical systems in the corresponding quantum case |l|J^|lO[ . 
Analyzing the Husimi projections of the eigenstates of the spin-boson Hamiltonian in relation to optical transitions we 
indicate a bridge between the eigenstate analysis of nonintegrable systems and observables in an optical experiment. 
In particular, we show how the intensity variations in the Husimi projections are related to the spectral randomness 
in the fine structure of absorption bands. Such spectral randomness in the eigenstates of the spin-boson Hamiltonian 
is also known as incipience of quantum chaos |id|-|l5|]. In this case random features of the spectrum are just appearing 
while the system is still far from displaying universal spectral fluctuations described by random matrix theory and 



known for excited states of spectra of polyatomic molecules |14|. 

A central point in our eigenstate analysis will be to find out to what extent the adiabatic reference systems of 
the spin-boson Hamiltonian are present in its exact eigenstates and how the appearance of spectral randomness can 
be interpreted as a mixing of such reference systems. Following jl3| we show that Bloch projections are a useful 
quantitative characteristic to describe this mixing. Computing the Bloch projections of the eigenstates, it is possible 
to distinguish the spectral region where the adiabatic branches of the spectrum are still intact from the region with a 
substantial mixing of adiabatic reference systems. Furthermore, by performing projections of the numerically obtained 
eigenstates onto the ground state we obtain the details of the fine structure of the absorption bands in the spectral 
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regions as indicated by the Franck-Condon principle. We note that the asymmetric case has some special features, 
which are not present in the symmetric situation: In the asymmetric case the final states of an optical transition can 
be located in the higher energetic regions of the corresponding adiabatic potentials. Then, as will be shown below for 
two representative cases, the fine structure of the spectrum and of the absorption bands can greatly vary from regular 
to irregular arrangements of absorption lines. 

The paper is organized as follows: In Sec. II the model and the basic equations are presented. In Sec. Ill the 
eigenstate analysis for two representative cases of the asymmetric spin-boson Hamiltonian is performed and in Sec. IV 
this analysis is connected with the properties of the absorption bands in the corresponding spectral regions. Finally, 
in Sec.V our conclusions are summarized. 



II. MODEL 



We consider a dimer system described by the Hamiltonian 

H = ^2(e n +lnQn)\n)(n\ + ^ V nm\n)(m\ + ^2 -(pi + ufal), (1) 



n n,m n 

n^rn 

where \n) (\m)) are the excited states of the molecular monomers constituting the dimer (n,m = 1,2) with the 
excitonic energies e n , the coupling constants to intramolecular vibrations 7„ and the transfer matrix element V nm . 
The variables q n , p n and u n represent coordinate, momentum and frequency of the vibrations, respectively. The 
coupling of the excited dimer states to an incident light wave with electrical field vector E n (t) is given by 

Hint = J2PnE n (t)\n}{0\, (2) 
n 

where jl n and E n (t) are the optical transition matrix element and electric field at the nth monomer, respectively, 
and |0) is the excitonic vacuum. We consider an asymmetric configuration of the dimer system by assuming different 
monomer energies e„ and optical transition matrix elements /2„. In order to keep the number of asymmetry parameters 
minimal we assume symmetric transfer rates, V nrn = V mn = —V, V > 0, and equal coupling constants and frequencies, 
i.e. 71 = 72 = 7 and oj\ = L02 = u). Then the coupling can be reduced to a single vibrational mode by introducing the 
relative displacement of the vibrations. 

The Hamiltonian (§) can be represented as an operator in the space of two-dimensional vectors C = (^) using the 
standard Pauli spin matrices (i — x, y, z). Passing to dimensionsless variables by measuring the energy in units of 
2V , H = H/2V, one obtains from ([!]) the spin-boson Hamiltonian 

H = e+t - ht x + i(P 2 + r 2 Q 2 )l + {J^rQ + e_)a z . (3) 

In (|^), the dimensionless relative displacement and corresponding momentum, to which the exciton is coupled, are 
given by Q — \/2V(qi — 92) and P = (pi — P2) /y2~V. The dimensionless parameters of the spin-boson Hamiltonian 
(p|) are related to the dimer system (Q) by 

r = W , (5) 

The parameter p expresses the coupling strength, and r is the adiabatic parameter measuring the relative strength 
of quantum effects of the subsystems. The sum e+ represents the center of the excitation energy, whereas the difference 
e_ is the asymmetry parameter. For e_ = 0, one obtains the symmetric spin-boson Hamiltonian with conservation of 
total parity (given by the operation Q — > —Q, <j z — > —a z ). For e_ ^ 0, this symmetry is broken. 

The adiabatic reference systems associated with (^|) are introduced by considering the eigenvalue problem of the 
adiabatic part /i a( j, 
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(7) 



The Q-dependent eigenvalues and eigenstates of h a & are easily obtained. In particular, one finds the two adiabatic 

ad ( 



potentials U^ d (Q) for the upper (+) and lower (-) states 




u>U(Q) = jQ \/t + ('-+ \/^) • !Si 

from which the Hamiltonians of the adiabatic reference systems are obtained 

H ± {Q) = \p 2 + U±{Q). (9) 

A key point in our analysis of the eigenstates and absorption bands will be the extent to which the exact eigen- 
states of the Hamiltonian (Q) can be viewed as a mixing of adiabatic reference systems connected with the adiabatic 
Hamiltonians As is shown below, the extent of this mixing can be controlled by analyzing Husimi projections 
and by computing Bloch projections from the numerically obtained eigenstates. Husimi projections constitute phase 
space representations of the eigenstates which can be compared with the phase space orbits of the adiabatic reference 
systems (|^). Bloch projections constitute another independent indicator of the mixing of the adiabatic branches in 
the spectrum. For the adiabatic branches these projections are calculated from the eigenstates of the Hamiltonian (]?]) 

\<ptd(Q)) = ^ (Vl + MQ)\ T> - y/l-MQ)\ I)) (10a) 
K d (Q)) = ^= (yi-A(Q)| T) + s/l + A{Q)\ [)) (10b) 



where 



Computing the expectation values x = (a x ) of the Pauli matrix a x for the adiabatic eigenstates (|10|), one obtains 
the Bloch projections of the adiabatic branches 



A(Q) = - 2e - + ^ . (11 ) 
y/l + (2e_ + VZprQ)* 



1 

y/l + (2e_ + V^PrQy 



(^ d (Q)\^d(Q)) = T 7t= . . (12) 



According to eq.(|l2|) there is a distinction in the sign of the Bloch projections for the adiabatic eigenstates. Com- 
puting the Bloch projections from the numerically obtained eigenstates this distinction in the sign will be used as an 
indicator for the case in which the spectrum resolves into adiabatic branches. 



III. EIGENSTATE ANALYSIS 

Our analysis is based on a diagonalization of the spin-boson Hamiltonian using a basis of product states of spin 
states with harmonic oscillator eigenstates. In this basis the eigenstates are obtained in the form 

w = E ( c ff it) + 4tl u>) m = E I s ™) h > ( 13 ) 

m m 

where |t) = Q, (||) = (?)) denote the spin up (down) states, \m) the harmonic oscillator eigenstates and cim (z =f, |) 
the expansion coefficients of the A eigenstate. For a fixed m we will also use the spin representation \ for the up 

and down coefficients cim (z =t,|) for a given oscillator state as indicated in the second part of ( |l3|) . The matrix 
dimension in the numerical diagonalization was N = 4000, for the eigenstate analysis the first 1100 eigenstates were 
used. 
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For the phase space representation of the eigenstates Husimi projections are used. Husimi projections h\ (a(Q, P); s) 
are defined by projecting an eigenstate |A) onto a set of coherent oscillator states \a(Q, P)), which scan the oscillator 
phase plane of the vibrational subsystem, while the spin projection \s) = o\ |f) + cj. |J.) is kept fixed, 

h x (a(Q,P);s)= |( A | «(Q, P), s)| 2 , (14) 



and \a(Q,P),s) = \a(Q, P))\s). Using the explicit form of |A) in (13) one finds 

h X [ni(). /-'):.si = ^( C S )c T + c £ )c l) Ha(Q,P)) 



where (m\a(Q,P)) = ^e"^ and a(Q.P) = ^(Q) + ^=(P>. 
For the Bloch projection of an eigenstate, 



(15) 



one finds from the expansion ( |l3| ) 



x x = {X\a x \\), (16) 



Before turning to the eigenstate analysis for definite parameter sets we note that the appearance of spectral ir- 
regularities can be expected in the high energetic region of the overlap of both adiabatic potentials only, where a 
sufficient number of states of both potentials mix. It is obvious that such a region can be probed by the final states of 
an optical transition in the asymmetric case only: In the symmetric case, as is evident from the adiabatic potentials 
(^) , the optical transitions would necessarily terminate in the ground state region of the upper potential, where only 
a few states of the upper potential can mix with the lower potential. Therefore when investigating the structure of 
the spectrum and eigenstates of the spin-boson Hamiltonian we paid special attention to a systematic change of the 
asymmetry parameter e_ combined with the coupling parameter p in order to produce appropriate configurations 
of the adiabatic potentials. Such configurations were found for relatively high coupling p and asymmetry values e_. 
Below we present the eigenstate analysis for two typical asymmetric configurations in the adiabatic parameter region 
(r < 1) with the parameter sets 

A: p = 4, r = 0.1 and e_ = 5 and 

B: p = 20, r = 0.1 and e_ = 10. 

The parameter set B corresponds to a stronger coupling value p and a larger asymmetry parameter e_, as compared 



to set A. To indicate the location of the absorption bands, which will be analyzed in section IV, the adiabatic potentials 
of set B and the ground state of the system are shown in Fig. [l]. 

In Fig. H we start with the eigenstate analysis for the parameter set A by displaying the Bloch projections x{E\) 
computed from the eigenstates according to cq. (|l7|) . In the low energy region one finds one adiabatic branch in the 
x{E\) dependence associated with the lower adiabatic potential (for which x\ > in accordance with the sign in (fL2|)). 
The overlap of the two potentials is marked by the appearance of a second adiabatic branch in the x(E\) dependence. 
The mixing of the eigenstates of the two adiabatic potentials is seen by the disappearance of both adiabatic branches 
and indicated by a broad band of Bloch projections in the high energy region. The presence of three characteristic 
regions in the Bloch projections with one adiabatic branch, two adiabatic branches and the mixing region is reflected 
in the Husimi projections of the eigenstates. In the regions where the adiabatic branches are intact one obtains regular 
Husimi projections concentrated around the phase space orbits of the integrable adiabatic reference Hamiltonians (|^) . 
Typical examples of such regular projections corresponding to equal spin projections c-f = cj, for the energetic overlap 
region with two adiabatic branches present are shown in Fig. ||(a) and Fig. |^(b). In Fig. @(a) the Husimi projection 
corresponding to an eigenstate in the upper adiabatic potential is shown, whereas in Fig. pTb) a projection associated 
with an eigenstates of the lower adiabatic potential is displayed. Sequences of such projections clearly attributable to 
the coexistence of two independent adiabatic branches were observed until the mixing region is reached. Then Husimi 
projections incorporating the phase space features of both adiabatic potentials are observed, as shown in Fig. |3](c). 

A closer inspection of these projections in the mixing region shows that although the distribution is concentrated on 
the regular phase space orbits of the adiabatic potentials, the intensity between these parts varies in an irregular and 
random way, when Husimi projections for a sequence of eigenstates are considered. The structure of the eigenstates 
which shows up in the random shift of intensity in the Husimi projections is a source of spectral randomness in the 
absorption bands. In the case of the parameter set A, however, this region of an irregular and random shift of the 
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phase space distribution is not reached by the absorption processes. In this case the absorption band is located in the 
region, where the two adiabatic branches are still intact and no spectral mixing between their eigenstates occurs. 

In the case of the parameter set B one finds a situation, where the spectral window, which is cut by the transition 
matrix elements of the absorption process, reaches the region of an irregular mixing of adiabatic states. The Husimi 
projections for selected eigenstates and equal spin projections which correspond to such a region and which are relevant 
as final states of an optical transition are shown in the set of Fig. [|(a)-(c). One observes a random variation of the 
intensity in the Husimi distributions of the eigenstates between the orbits of the two adiabatic potentials. As we show 
in the next section, as a consequence of this random variation intensities of lines for optical transitions terminating 
in these states vary randomly. 



IV. ABSORPTION 



In the calculation of the optical absorption the dimer system is assumed to be small compared to the optical 
wavelength. Then the electric field vector is approximately equal at the monomer sites of the dimer configuration, 
i.e. E n (t) = E(t). Introducing the projections //„ of the monomer transition matrix elements onto the direction of 
the field polarization, i.e. E(t)fl n — E(t)[i n , (n = 1,2), the interaction of the dimer with the incident light wave is 
represented in the form 



Him = E(t)M, 



int •> 



(18) 



where Mi nt is the interaction operator 



= ^\n)n n ( 



(19) 



In the ground state there is no excitonic-vibronic interaction, i.e. the ground state of the system \g) is given by the 
direct product of the excitonic vacuum |0) = |ea;cO) and the vibrational ground state. For the ground state of the 
vibrational subsystem we assume the zero temperature case, where it is in its lowest state |m = 0), with |m = 0) the 
rn = Hermitian polynomial of the undisplaced Q oscillator. Then \g) = \exc0)\m — 0). For the optical transition 
matrix element M\ g = (\\Mmt\g) of the interaction operator (|l9|), which corresponds to an absorption process from 
the ground state g into an excited state A, one then obtains using eq.(13) the expression 



M 



As 



(A) 

/xic T0 



(A) 



(20) 



The matrix element (BG) is representable in the form of a particular projection of the spin vector associated with 



the m = vibrational state in (|l3|), 



4 A) ). 



Introducing the angles 



cos a 



Ml 



i4 



, sin a 



M2 



(21) 



and a spin vector \su) defined by 



cos a 
sin a 



(22) 



one represents 



in the form 



(23) 



Measuring the square of (|23|), i.e. the absorption strength of an optical transition from \g) to |A), in units of 
(fii + /tzi) one finds for the dimcnsionsless absorption strength Q\ g 



M l 9 



(24) 



Comparing eq.(|24|) with the expression for the Husimi projection (15) it follows that Q\ g is equal to a Husimi 
projection of the A-eigenstate of (||) taken at the phase space point Q = P = with the spin projection being fixed 
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at : For Q = P = one obtains a(Q, P) = 0, which selects the m — term in the sum of (|l5|), inserting the value 
s M for the spin vector s, one finds that /ja(0; s) reduces to the r.h.s. of eq. (|24|) , i.e. 

Qxg = h x (0;s fi ). (25) 

Eq. ( p5[ ) is in line with the Franck-Condon principle for optical transitions: The transition is vertical from the 
ground state region at Q = P = in the phase space representation and probes the intensity of the final state by 
its Husimi distribution in the same region. The optical transition matrix elements of the monomers define the spin 
projection s M . For the special case of a symmetric dimer with /xi = [12 the transition occurs into the symmetric 
combination of monomer states. Differences in the optical transition matrix elements, i.e. [i\ ^ ^2, introduce a 
second asymmetry parameter besides that of the site energy asymmetry e_ . These differences in the optical transition 
matrix elements express the optical asymmetry of the dimer and enter the optical dimer matrix element ( pof ) through 
the spin projections j2l|). 

In Fig. ||(a) the absorption bands Q\ g calculated from the numerically obtained eigenstates for the parameter set 
A are represented as stick spectra. Absorption bands are located in accordance with the Franck-Condon principle. 
This is seen by comparing their positions with the final states following from the adiabatic potentials in <Js|) by setting 
Q = 0. One finds an energetically lower band for optical transitions terminating in the lower adiabatic potential and a 
higher band for transitions occuring into the region of the overlap of the two potentials. The shape of the lower band 
on the left side of Fig. |^(a) is completely regular. The lines of the higher band on the right side of Fig. ||(a) are also 
regularly arranged. A closer inspection of the higher band, however, reveals that between the lines shown in Fig. |^(a) 
lines with a much weaker intensity are embedded, which form another regular band. This is evident from the inset, 
Fig. ||(b), where parts of the lines of the weak and strong bands are shown together. In order to make the weak band 
visible in the inset a much smaller intensity scale is used. The two bands, into which the absorption spectrum of 
the high energy part resolves, are easily identified as to belong to the states of the two coexisting adiabatic branches 
analyzed in the preceding section for the parameter set A: The weak band with relatively small intensities is due to 
optical transitions terminating in the high energy states of the lower potential (the sign of the Bloch variables of the 
final states is x\ > in accordance with the sign of the lower adiabatic branch ([l2|)), whereas the strong band with 
much greater intensities is due to optical transitions into the low energy states of the upper potential (for these states 
x\ < 0, again in accordance with the sign of the adiabatic branch (|12j)). 

In Fig. ^| the absorption bands are shown for the parameter set B. The bands are located in the spectral regions 
of the Franck-Condon energies indicated by the arrows in Fig. [I]. The lower absorption band, shown on the left side 
of the Fig. §, is regular like the lower band in the Fig. ||(a). The higher band, displayed on the right side of Fig. ||, 
however, is completely irregular and cannot be resolved into independent subbands as in the case of the parameter set 
A. This is a consequence of the mixing of the adiabatic reference systems and the irregular structure of the eigenstates 
analyzed by the Husimi projections in the preceding section. In particular, these projections (see Fig. 0) show a 
random variation of the intensity in the Q = P = region, which according to eq.(p5|) is relevant for the final states 
of the optical transitions. This random variation is probed by optical transition matrix elements and results in an 
irregular pattern of lines in the fine structure of the high energy absorption band. 

The absorption bands in Fig. || and Fig. [6] were calculated for equal optical transition matrix elements at the 
monomers fi± — (12, which correspond to spin projections of the eigenstates with equal components. The changes 
occuring in the lower and upper absorption bands in the limiting cases of an optical asymmetry with n± ^ 0, ^2 — 
and [i\ =0,^2^0 for the case of the parameter set B are compared in the Figs. [?] (a) and (b), respectively. These 
limiting cases correspond to an optical asymmetry, when one of the monomers constituting the dimer is optically 
active only. Then according to eq.(|2^) the optical transition matrix element is determined by either the spin up ci^ 

or the spin down c9q coefficients of the A eigenstate. As is seen from Fig. ^|, which is representative for the case of a 
positive sign of the excitation energy asymmetry e_ (e_ = 10 > in the case shown in Fig. ^), the lower and upper 
bands show a redistribution of their intensities with the optical asymmetry: For the given sign of the excitation energy 
asymmetry e_ > (e_ < 0) the intensity of the upper band is increased (decreased) compared to the symmetric case, 
Mi = displayed in Fig. k| The lower band behaves in the opposite way, the peak intensities of both bands differing 
by a factor of about 10 -3 . In particular, maximum absorption is reached for the high energy band, when the monomer 
with the higher excitation energy is optically active only, i.e. e_ > with fj,i ^ and [12 = 0. Reversing the optical 
activity of the monomers but still considering the same asymmetry in the excitation energies, i.e. e_ > with \i\ = 
and y,2 7^ 0, one finds maximum absorption for the low energy band (if e_ < 0, the quantities [i\ and /i2 have to be 
interchanged in the above consideration). As in the case of equal transition matrix elements the lower bands shown 
on the left sides of Fig. |^(a) and Fig. f7](b) display a regular fine structure, whereas the fine structure of the upper 
bands, displayed on the right sides of Fig. ^(a) and Fig. f?](b) is irregular. We note that the appearance of the weak 
bands is due to the coupling of the excited states of both monomers: The monomer, which is optically not active, is 
still present in the absorption spectrum due to the coupling of the excited states in the dimer system. In the case 
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displayed on the right side of Fig. 0(b), e.g. an optically nonactive monomer with a high excitation energy attached to 
an optically active monomer with a lower excitation energy produces a weak but irregular band in the dimer system. 
A closer inspection of the fine structure of the irregular bands for the limiting cases of an optical asymmetry, in 
Fig. |(a) and Fig. |(b), together with the irregular band for the symmetric case, shown in Fig. ||, reveals that the 
sequence of lines with a strong and small absorption strength is almost identical in all bands, there is only a small 
overall change in the line intensities in all three representations. The reason behind this behavior becomes evident 
after an inspection of the ratio of the spin projections 

rW (26) 

which is presented in Fig. ^| as a function of the eigenstate energy E\ for the spectral regions corresponding to the 
lower and upper absorption bands, respectively. As is evident from Fig. |§| is a smooth function of the energy 
E\ in the spectral regions of both bands. In particular, for the case of the irregular upper absorption bands, for 
which in a sequence of eigenstates the spin projection coefficients and cIq vary in an irregular way, their ratio 
r( A ) considered between neighbouring eigenstates is practically identical. This correlation between the spin up and 
spin down (c|.q and c^, ) coefficients results in an almost smooth interpolation of the fine structure of the absorption 

bands between the limiting cases displayed in the Fig. The smoothness of the function can be used to express 
the absorption strength for arbitrary optical asymmetry by one of the limiting cases. For example, representing the 
down coefficients c[q by the up coefficients ci^ in (^o| ) using eq. (|26l), one obtains 

PI "T / i 2 

where Q\ g = [c^] 2 is the absorption strength for the case /ii ^ and /12 = 0. The smoothness of the ratio of the 
different spin projections together with eq. (|27]) allows an interpolation between absorption bands for different 
optical asymmetries. In particular, using this ratio in the prefactor in front of Q~^ g the relative intensities of absorption 
bands corresponding to different optical asymmetries can be determined. 



V. CONCLUSIONS 



The transition from the coexistence to a mixing of the adiabatic branches in the eigenstates of a symmetry-broken 
spin-boson Hamiltonian can be controlled by Bloch and Husimi projections and it shows up in the structure of 
absorption bands in excitonic-vibronic coupled dimer systems with asymmetric adiabatic potential configurations. 
In the mixing region of the adiabatic branches phase space distributions of the eigenstates are concentrated on the 
phase space orbits of both adiabatic potentials and their intensity varies in a random and irregular way. This random 
variation in the intensitiy distribution of the eigenstates is probed by the optical transition matrix elements and 
it is a source of spectral randomness in the fine structure of the absorption bands. The asymmetry of the optical 
transition matrix elements of the dimer configuration is representable by suitably chosen spin projections of the spin 
part of the eigenstates. In the limiting cases of the optical asymmetry, in which one of the monomers constituting 
the dimer is optically active only, the intensities of the absorption lines are determined by either the spin up or spin 
down coefficients. For the states relevant for the optical absorption correlated spin down and spin up coefficients are 
obtained. As a result the shape of the fine structure of the irregular absorption bands remains almost the same for 
different optical asymmetries. Finally we point out that in the case of a symmetry-broken spin-boson Hamiltonian 
optical absorption processes can probe the high energy regions of the overlap of adiabatic potentials where it is possible 
to observe a transition from regular to irregular absorption spectra for asymmetric molecular dimer configurations. 
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FIG. 1. Asymmetric adiabatic potentials and ground state configuration for the parameter set B. In the case of the parameter 
set A the asymmetry in the configuration of the adiabatic potentials is similar, but less pronounced as compared to B. Arrows 
mark the location of the absorption bands at the Franck-Condon energies. 



8 



0.3 




FIG. 2. Bloch projections x{E\) of the eigenstates for the parameter set A. The presence of three characteristic regions with 
one adiabatic branch, coexistence of two adiabatic branches and the mixing region is clearly visible. The Bloch projections for 
the parameter set B show a similar behavior. 
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(a) 



(b) 



(c) 




FIG. 3. Husimi projections for the parameter set A in the overlap region of the adiabatic potentials. The solid lines indicate 
the classical phase space orbits of the adiabatic Hamiltonians (^). The projections correspond to the eigenstates 120 (a), 121 
(b) and 310 (c), with the eigenstate energies E\20 = 5.00 28, E\2\ = 5.0956 and £310 = 14.5176, respectively. The selected phase 
space parts are (a): -50 < Q < 50, -5 < P < 5, (b): -75 < Q < 75, -6 < P < 6 and (c): -90 < Q < 90, -8 < P < 8, with 
the Q-axis displayed horizontally and P-axis displayed vertically. The projections in (a) and (b) are located on the phase space 
orbits of the upper and lower adiabatic potentials, respectively, and show the coexistence of both adiabatic branches in the 
selected energy interval. The projection (c) is located on the phase space orbits of both adiabatic potentials and characteristic 
for the mixing region of the adiabatic reference states in the spectrum. 
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(a) 




FIG. 4. Husimi projections for the parameter set B. The solid lines indicate the classical phase space orbits of the adiabatic 
Hamiltonians (^|). The projections correspond to the eigenstates 300 (a), 302 (b) and 303 (c), with the eigenstate energies 
-E300 = 10.0237, -E302 = 10.1467 and -E303 = 10.1520, respectively. In all projections the selected phase space parts are 
— 120 < Q < 120, — 10 < P < 10, with the Q-axis displayed horizontally and P-axis displayed vertically. The projections are 
located on the phase space orbits of both adiabatic potentials and characteristic for the mixing of adiabatic reference states. 
One observes a random variation in the intensity of the distribution between the different eigenstates. This random variation 
includes the center of the rectangle, i.e. the region around Q = P = 0, which is the final state region for the upper absorption 
band shown in FigM The different intensities in this region correspond to different strengths of the absorption lines of the 
upper band in Fig.HT 
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FIG. 5. Lower and upper absorption bands for the parameter set A. In (a) the strong intensity lines of the upper band are 
visible only. A closer inspection of the upper band shows that this band is a superposition of two regular bands corresponding 
to the final states in the upper and lower potential and the coexistence of two independent adiabatic branches in this part of 
the spectrum. This superposition becomes evident from the inset (b), in which a change of scale is used to display the weak 
lines, which are embedded between the strong lines (not drawn to peak intensity and cut off at the upper edge). In the inset a 
pair of neighbouring weak and strong lines are indicated by the labels w and s, respectively. 
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FIG. 6. Lower and upper absorption bands for the parameter set B and the case of optical symmetry, pi = ^2- A broken 
energy scale is used to display both bands. The lower band is regular with final states in the lower adiabatic potential, whereas 
the upper band is irregular due to the mixing of the adiabatic reference states, compare with the intensity variation in the 
Husimi projections displayed in Fig. 0. 
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FIG. 7. Absorption bands for the parameter set B in the limiting cases of optical asymmetry, when one of the molecules 
of the dimer is optically active only: Hi ^ 0,1*2 = (shown in part (a)) and /t*2 0, fii — (shown in part (b)). Note the 
differences between the overall intensities of the lower and upper bands as indicated by the scales on the left and right hand 
sides, respectively: In (a) the intensity of the lower band is by three orders of magnitude smaller than the upper band, in 
(b) the intensities of the bands are reversed. Independent of this change in intensity the upper bands in both (a) and (b) are 
irregular and have a similar fine structure. 
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FIG. 8. Ratio of the spin down to spin up coefficients r'"* 1 ' in the spectral regions of the lower and upper absorption bands 
of Fig. fij. The scales on the left and right side correspond to the ratios in the region of the lower and upper band, respectively. 
Note the smooth dependence of r y ' on the eigenstate energy in the spectral region of the upper absorption bands, despite the 
irregular structure of the upper bands in Fig. Hand Fig. M (a),(b). 



15 



